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The evolution of scale-invariant gravity waves from the early universe is analyzed 
using an equation of state which smoothly interpolates between the radiation domi- 
nated era and the present matter dominated era. We find that for large wavenumbers 
the standard scale-invarant wavefunction for the gravity wave severely underestimates 
the actual size of the gravity wave. Moreover, there is a definite shift in the phase 
of the gravity wave as it crosses the radiation-matter phase transition. The tensor- 
induced anisotropy of the cosmic microwave background and the present spectral 
energy density of the gravity wave is then calculated using these results. 
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It is well known that a scale-invariant (SI) spectrum for gravity waves (GW) can be 
produced in an inflationary cosmology UJ. During inflation, fluctuations in the gravitational 
field will be red-shifted out of the horizon, after which they freeze and remain at a constant 
amplitude. Much later when they re-enter the horizon during either the radiation-dominated 
(RD) or matter-dominated (MD) era, they will appear as classical GW. Via the Sachs- Wolfe 
effect H, these primordial GW will leave imprint on the cosmic microwave background 
(CMB) ||. In fact, recent measurements of the power spectrum extracted from the large- 
angular-scale anisotropy of CMB [[| suggests a behavior that may be explained by these 
primordial GW ]5j|| . Moreover, the direct detection of GW is one of the main goal of many 
future experiments. Tracing the precise evolution of these GW is therefore mandatory and 
important. 

In this Letter we shall analyse in detail the time evolution of a SI spectrum of GW in 
the standard hot big-bang cosmology. This subject has been extensively studied before, but 
in all previous approaches one either considers only the long wavelength behavior of the 
GW, and/or a sudden change of equation of state as the universe moves from the RD to 
MD eras is adopted. Here the novelty is that we shall include a smooth radiation-matter 
phase transition, which will allow us to obtain specific solutions of the wave equation across 
the transition. By doing so, we shall find a spectrum for the short-wavelength GW that is 
radically different from what is currently accepted. Moreover, we shall be able to explicitly 
determine the phase of the GW. We then apply our result in the calculation of the tensor- 
induced anisotropy of CMB, and the present spectral energy density of GW. Attempts at 
incorporating a smooth RD to MD transition in the equation of state has also been made 
by using transfer function methods |7[]. 

We assume a spatially flat two-component universe containing radiation and dust. 
The metric that we will use throughout is the flat Robertson- Walker metric: ds 2 = 
dt 2 — a 2 (t)dx. 2 = a 2 (r]') (dr]' 2 — dx 2 ) where a(rj') and drj' = dt/a(t) are the scale factor and 
conformal time respectively. We will consider a monochromatic plane GW of wavenumber 
k' . We first define the new variables t] = (\/2 — l)rj' /rf eq , k = k'r]' eq /(\/2 — 1). rj' eq is the time 
at which the energy density of radiation is equal to that of matter. The equation of motion 
of the GW amplitude is then given by || 
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In a generic inflationary model, a SI spectrum of GW can cover a very wide range of k (for 
instance, 10~ 5 < k < 10 23 ) ||. For these fc's the standard boundary conditions for h{rj) are 

h(rj) = 1 and = as 77 -> 0. (2) 

a?7 

In the standard solution of a SI spectrum, which we shall call the SI solutions, Eq. ([!]) 
is solved in the RD and MD universes seperately 0] giving 

, , v { 3o(kv), if V^Veg] 

hrmiV) = < n . (3) 

where j u {kr]) and y u {kr]) (below) are the spherical Bessel functions of order i/ of the first 
and second kind, respectively. 

As a first attempt to take into account the radiation-matter phase transition in the 
analysis of the evolution of G W, one typically makes the sudden approximation 

2 da ( 2 /v, if V < Veg] 
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{ 4/(77 + Veq), Hv>Ve g . 
Once again it is straightforward to solve Eq. ([]]) giving 

(joM, if 77 < 77 eg ; 

= < (5) 
I (Veq/V) (AsMkn) + B s yi(kr])) , if 77 > r/ e? . 

The subscript s denotes the fact that we are using the sudden approximation (||) . Notice 
that h s (j]) satisfies the boundary conditions (0), while the integration constants A s and B s 
are determined by requiring h{rj) and dh/drj be continuous at 77 = 77^. This gives 



cos krj eq sin kr\ eq sin 3 A;77 e 
2ki] eq 4:(kr] eq ) 2 2(krj eq) 



( sin kr] eq cos kr] eq cos 3 kr] eq \ 
B s = -4 1 cos kr] eq - —- — - + -| . (6) 



leg) 



A s -A[ sin /••//,„ + + , 2 . ^ 2 

^3 

2%, 2(A;77 e9 ) 2 ' 2(k Vt 
Notice also that h s —>■ 1 as krj — > 0, and by this criterion /i s is a scale-invariant solution. 
Similar solutions can also be found in Ref. JTH ]. 

Of course, the evolution of the universe is not nearly this sudden. There is instead a 
gradual change in the equation of state from the early RD era into the MD universe of the 
present day with the scale factor 



a{v) = a eg V (v + 2) • ( 7 ) 

a eq is the scale factor at rj = i] eq . (Note that r\ eq ~ 0.41 and the present time rjo m 76.5 if we 

take a /a eq = 6000.) Then Eq. (|l|) reduces to 

^ + 4(,+ l)rift + ^ = 
drj 2 r](r) + 2) ar\ 



V = o, (10) 



This equation can be solved exactly in terms of spheroidal wave functions |ll|jr2|] . Doing so 
will not be very illuminating and we shall instead use the WKB approximation to obtain an 
approximate analytical expression. 
First, we define 

h(v) = gfo) ( 9 ) 

Then Eq. (f|) becomes 

2_ 

dr\ 2 \ 77(77 + 2) 

which is analogous to the Schroedinger equation for the "wavef unction" of the particle with 
energy k 2 in a "radial" potential ^(77) = 2/{i](ri + 2)} since rj > 0. Because the boundary 
conditions for h(rj) are given at 77 = 0, the problem now reduces to the quantum mechanical 
problem of a particle with energy k 2 tunneling out of the potential V{rf). 

Notice that in this language the sudden approximation corresponds to a particle tunneling 
out of the potential 

f 0, if 77 < r] eq: 
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The singularity in V(rj) at rj = has been removed and has been replaced by a potential 
well with width 7] eq and height 1/2. This would explain the resonances in the amplitude of 
h s (rj) when kr) eq ~ nn for some integer n, as can be seen explicitly in Eq. (§). 

Since 77 > 0, to apply the WKB approximation define 77 = e x where x G (—00, 00). Then, 
taking y(x) = e x ^ 2 u(x), 

— + e 2a; r - - u = 0. 12 
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From the standard WKB connection formulas across the turning point, we find that for 

7/ < 77t, 
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while for 77 > rjT, 
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The tunneling amplitude A(k) is 
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and rjT is the turning point of the potential in Eq. (0) which is determined by 
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One can show that 1/2 < /c?7t < 3/2. We should also note that when kt] is large, and when 
i] ^> 2, the transfer function T(fc) of Ref. can be easily obtained from the WKB solution: 
T{k) w (2/3)A; 3 / 2 A(A;). 

We shall now compare these three approximate solutions with each other. Taking the 
small kr] << 1 limit corresponding to modes well outside the horizon, we find that h rm (r)) ~ 
h s (rj) ~ 1 and these solutions are scale invariant. For the WKB solution, on the other hand, 
we find that h w {rj) w H(r]), an 77 dependent function that is close to unity for all 77. This 
unexpected //-dependence is mainly due to the breakdown of WKB approximation in this 
limit. Since h satisfies the boundary conditions (|2]), despite this small deviation from unity 
the wave amplitude will remain constant when it is outside the horizon. 

Next, the asymptotic large 77 ^> 1 limit for k fixed at some finite value gives 
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for the MD equation of state while for the sudden approximation, 

-3 cos kf]/(kr]) 2 if k 1; 
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The WKB approximation gives 



^0- cos(A;?7 - 7r) /(krj) 2 if k < 1; 
^(77) « <; v ( 19 ) 
2 x Acos(fcr/ - 7r/4)/(/b/ 2 ) if Jfe > 1. 



For small fc all three results reduce to Eq. fll7l) , except a slightly different amplitude and an 
extra phase of it in This is to be expected. When the wavelength of the GW is much 
greater than r] eq , it does not "see" the RD era and is affected only by the MD universe. 
Consequently, the SI solution h rm in the MD era works very well for GW of wavelength in 
the regime kr] eg <C 1 at late times. 

The situation changes dramatically for short-wavelength gravity waves, however. Denot- 
ing the amplitudes of h rm (r]), h s (rj) and h w {rj) by M rm , M s and M w respectively, 
M rm 3 M rm 3 fe M s „ fe 



M s 4kr] eq ' M w 2k\ 2 ' M w 2r]eq \/ 2 ^ 
for k ^> 1/4. Since the WKB solution is a closer approximation of the actual history 
of the universe, we see that the sudden approximation is a very good approximation of 
the amplitude of the GW. The SI solution in the MD era, on the other hand, severely 
underestimates the amplitude of the GW when k > 1.75. 

Although the change in the amplitude is most drastic for large k, the relative change in 
the phase of h{rj) between the SI solution and the sudden and WKB approximations are 
present at all k and may be significant fl3|]. In fact, we see that although the amplitude of h s 



very closely approximates that of h w , its phase is k dependent which, more importantly, does 
not approach any fixed value for large k. Thus, for different values of k the phase difference 
between h s and h w will differ. These phase differences can be seen explicitly by comparing 



Eqs. (|T^-([nj) and also in Fig. 1. Here we have also plotted the numerical solution h n (r]) 
to Eq. (§. 

In our sudden approximation we could have also approximated the equation of state with 
an a which has a discontinuous first derivative; namely da/dr] = 2a/r] for 77 > r] eq . In this 
case we would have found an h s which, for k 3> 4, underestimates the amplitude of the GW 
by a factor of 4. Its phase difference, however, approaches a constant value of 7r/4 at large 
k. 

We shall now use these results to calculate the anisotropy of CMB induced by the SI 
spectrum of GW as well as its spectral energy density. From the Sachs- Wolfe effect, the 



formula for the power spectrum G% is given in, for example, ||. Integrating from the decou- 
pling time r] dec ps 1.54 for a /arf ec = 1100 and taking the inflation parameter v = Vo/raj^, 
where Vo and rripi is the de Sitter vacuum energy and Planck mass respectively, we have 
numerically evaluated C\ by using both the SI solution h rm and the numerical solution h n . 
A comparsion of the two results is given in Table 1. As expected, for small I the differences 
are insignificant since at these values the dominant contribution to C\ come from small k 
modes, precisely where the SI solutions work well. Due to the difference in the factor of 1/k 
in the amplitudes of the SI and WKB solutions there is a change in C\ for large I, but unlike 
the results from using transfer functions J7|, we find that the SI solution overestimates the 
numerical solution around / ~ 100. This, we believe, is due to the phase shift between h n 
and h rm which cannot be taken into account by using a simple transfer function. 
As for the spectral energy density of the GW 

a,= £ (2D 

A=+,x Pc aK 

which is a measure of the total amount of energy deposited into the experiment at any 
frequency k. p c is the closure density of the universe. To estimate the amplitude of this 
quantity, one typically makes use of the following physical argument. Before the GW enters 
the horizon, it is a pure scale- invariant wave with h = 1. After it enters the horizon, on 
the other hand, the GW behaves effectively as radiation and will scale with a as such. This 
gives a tt g ~ 10" 13 for v ~ 10" 9 0,|. 



We can also calculate fl g directly, however. Since 

£^3^J%>| 2 , (22) 

where G and fc p h ys are the gravitational constant and physical wavenumber respectively, once 
h{rj) is known, one can simply evaluate fi g (r]) at the present time rj . In fact, using Eq. ([5J, 
we find that 

^M = ^>(r/o)| 2 f^r) 2 - ( 23 ) 

97r y T] + 1 ) 

Since t]q 3> 1 and since any experimentally detectable k is extremely large in our units, we 
find that by using the SI solutions in Eq. (p3|), 
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V g (Vo) w = cos2 (^o - vr/4), (25) 



2v 1 

ttJri ) = — - — cos 2 (kri ). (24) 

Because k ~ 10 9 , the amplitude of Q g (r]o) rm is infinitesimally small, contradicting the above 
physical argument. If, on the other hand, we use either sudden approximation solution or 
the WKB solution, then 

32v fri \ 2 

^a(Vo) s = ^ cos 2 (kr] -n/2-kr] eq ), 
Wve^ 1 1 

97T 7?g 

and the amplitude of both Q g is a constant independent of k. 

Since Q g (r)) is an oscillatory function, we have plotted in Fig. 2 the time average 

k rvo 

n g = — n g ( v )d V} (26) 

Z7T Jj7o-27r/fc 

of fi 9 over a complete period at the present time verses k for k much larger than the horizon 
size ku = 27r/?7o using the four solutions of (|l|). It is, moreover, this time average which is 
measured in experiment and not Q g {i]). For k near kfj, on the other hand, we have plotted 
Q g (rjo) itself evaluated at the present time 770 verses k in Fig. 3 since for these values of k, 
Q g oscillates too slowly for the averaging to be meaningful. Once again we can see from 
Fig. 2 that the fi 9rm obtained using h rm severely underestimates the actual Q 9n obtained 
using h n for large k. As expected, Q 9s and Q 9w obtained using h s and h w are both fairly 
good approximations of VL 9n . Notice also the dip in fi 9 (?7o) in Fig. 3 when k ~ 0.7kn- This 
sharp decrease in Q 9 at the present horizon size is due the finite size of the universe. 

In summary, we have studied the temporal dispersion of a SI spectrum of primordial GW. 
We have found that the standard SI solution used in the literature severely underestimates 
the size of the short wavelength, scale-invariant GW at the present time. Although we 
would thus expect the C/ calculated from the SI solution to always underestimate the actual 
C\ calculated numerically, this does not turn out to be the case. For / ~ 100, it will 
overestimate the actual C\ and is due to the relative phase shift in the h n verses h rm . This 
underestimation by the SI solution also occurs with the calculation of Q 9 , which up to now 
has been calculated using either indirect physical arguments or by using transfer function 
methods 0. This problem does not exist when using either the sudden approximation or 
the WKB solutions, however, and physically reasonable results for Q g can be obtained. 
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We also see from our solutions that the GW picks up a shift in phase as the universe 
undergoes the radiation to matter phase transition which is completely overlooked not only 
in the standard SI solutions, but also in the transfer function method. This phase shift of 
the GW may conceivably be an invaluable probe of the early universe before the hydrogen 
recombination. We can imagine future experiments performing interferometry on two inci- 
dent GW coming from different parts of the sky with their phase shifts being dependent not 
only on the coincident wavenumber k, but also on the detailed structure of the radiation- 
matter phase transition. More importantly, it can also serve as a probe of the anisotropy 
of the early universe as, conceivably, one part of the universe could have undergone the 
phase transition before another part, thereby producing a relative phase shift. However, 
the presence of initial randomly distributed phases will complicate this situation, as they 
should appear in the form of random noise in the measurement of the phase shift. Such 
experiments would require much greater sensitivity or new techniques than is now possible, 
of course, and would lie in the far future. 
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Captions 



Table 1. Table of C\ values calculated with the standard SI solutions, denoted by Ci rm , 
and the numerical solution of Eq. denoted by C%. 

Figure 1. Graphs of h rm (r]), h s (rj), h w {jj) and h n {rf) to Eq. (||) for k = 200k H where 
k H = 27T / 770 . Notice that among the three approximate solutions, only the WKB solution 
accurately reproduces the phase of h n . 

Figure 2. Graphs of Q grm , Q gs , Q gw as well as Q 9n verses k for large k. Notice the drastic 
drop in Q n 

Figure 3. Graphs of Q g (r] ) rm , fl g (r] ) s , Q g (r] ) w as well as fl g (r] ) n verses k for small k. All 
three approximate analytical solutions are good approximations of h n in this regime. The 
dip in Q g (r]o) at k « 0.7kn is actually a point at which Q g vanishes and is due to the current 
finite size of the universe. 
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Table 1 
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